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We address the following question: Which kinds of symmetry protected topological (SPT) Hamil¬ 
tonians can be many-body localized? That is, which Hamiltonians with an SPT ground state have 
finite energy density excited states which are all localized by disorder? Based on the observation 
that a hnite energy density state, if localized, can be viewed as the ground state of a local Hamil¬ 
tonian, we propose a simple (though possibly incomplete) rule for many-body localization of SPT 
Hamiltonians: If the ground state and top state (highest energy state) belong to the same SPT 
phase, then it is possible to localize all the finite energy density states; If the ground and top state 
belong to different SPT phases, then most likely there are some finite energy density states which 
can not be fully localized. We will give concrete examples of both scenarios. In some of these 
examples, we argue that interaction can actually “assisC localization of finite energy density states, 
which is counter-intuitive to what is usually expected. 


1. INTRODUCTION 

Symmetry protected topological (SPT) states and 
many-body localization (MBL) are two striking phenom¬ 
ena of quantum many-body physics. A d—dimensional 
SPT state is the ground state of a local Hamiltonian 
whose d—dim bulk is fully gapped and nondegenerate, 
while its (d — 1)—dim boundary is gapless or degenerate 
when and only when the system preserves a certain sym¬ 
metry G mm- An SPT state must have “short range en¬ 
tanglement” ; meaning that the entanglement entropy of 
its subsystems scales strictly with the area of the bound¬ 
ary of the subsystem: Sa ^ L'^~^ [5], where L is the 
linear size of the subsystem A. MBL refers to a phe¬ 
nomenon of the entire spectrum of a local Hamiltonian 
with disorder, including all of the highly excited states 
with finite energy density. Localization of single parti¬ 
cle states under quenched disorder is well-understood [1] , 
and recent studies suggest that localization can survive 
under interaction 011. In our current work the phrase 
MBL refers to systems whose all many-body eigenstates 
are localized, namely the entanglement entropy of all fi¬ 
nite energy density states obey the same area law as SPT 
states instead of the usual volume law typically obeyed 
by finite energy density states. 

These observations imply that in a many-body local¬ 
ized system, any finite energy density state actually be¬ 
haves like the ground state of a local parent Hamiltonian. 
Indeed, it was proposed that phenomena such as stable 
edge states and spontaneous symmetry breaking 00- 
0, which usually occur at the ground state of a sys¬ 
tem, can actually occur in finite energy density states 
of MBL systems. In fact, we can define a MBL sys¬ 
tem as a system for which any finite energy den¬ 
sity eigenstate is a short range entangled ground 
state of a local parent Hamiltonian. And if the 
system preserves a certain symmetry, then any finite 
energy density state of the MBL system should 
also obey the classification of SPT states. Then we 
can view energy density e as a tuning parameter between 


SPT states. Of course, in the thermodynamic limit, be¬ 
cause there are infinite states in an infinitesimal energy 
density interval (e,e -|- de), we expect there exists many 
Id curves in the spectrnm parameterized by s with one 
state |'0)e at each e, which is the ground state of an ef¬ 
fective SPT Hamiltonian H^. And on each such curve 
|'/)e is (roughly speaking) continuous in the sense that 
|'/)e and \ip)e+de are similar (despite being orthogonal), 
namely physical quantities averaged over the entire sys¬ 
tem change continuously with e on this curve. In an 
ergodic system, the eigenstate thermalization hypothe¬ 
sis m implies that most states with similar energy den¬ 
sity e are similar (their reduced density matrices all be¬ 
have like a thermal density matrix); in a MBL state, 
although states with the same energy density can in prin¬ 
ciple be very different, we still expect (assume) that the 
continuous curves mentioned above exist, although states 
in different curves can be very different. 

Within one of these curves mentioned above, tuning 
e is just like tuning between the ground states of local 
Hamiltonians. Furthermore, by tuning e there may or 
may not be a phase transition. In particular, if all excited 
states belong to the same SPT phase for arbitrary energy 
density e, then there does not have to be any quantum 
phase transition when tuning e, which implies that all 
of the excited states have short range correlations and 
area-law entanglement entropy, i.e. all the finite energy 
density states are localized; on the other hand, if states 
with different energy density e on the same curve belong 
to different SPT phases, then there must be at least one 
phase transition at certain critical energy on this curve 
when tuning e. This phase transition behaves just like 
an ordinary zero temperature quantum phase transition 
between different quantum ground states under disorder. 
For Id systems this “critical” energy density state could 
be in the “infinite-randomness” phase EHH], whose en¬ 
tanglement entropy scales logarithmically with the sub¬ 
system size m, hence it is not fully localized. The exis¬ 
tence of the “infinite-randomness” states at finite energy 
density have already been observed in Ref. 1161 
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Due to the fact that in a generic nonintegrable Hamil¬ 
tonian H, the ground state |G) and top state |T) (highest 
energy state of H and also ground state of —H) are usu¬ 
ally the easiest states to analyze, the most convenient 
way to determine the existence of “critical” states in the 
spectrum is to check whether the ground and top states 
belong to the same SPT phase or not. In summary, if |G) 
and |T) belong to different SPT phases, and if we under¬ 
stand that these two SPT states are separated by one or 
multiple continuous phase transitions (this will depend 
on the type of SPT phases |G) and \T) belong to), then 
there must be some “critical” excited states in the spec¬ 
trum which cannot be fully localized |46j . We will apply 
this rule to various examples in the next section. 

2. EXAMPLES 

2A. Kitaev’s chain: localization 

We first apply our argument to the Kitaev’s chain: 

H = + {-ySt + Atj) ijjjj+i, (1) 
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where are Majorana fermions and Atj is a random 
hopping parameter with zero mean and standard devia¬ 
tion CTAt- The topological superconductor phase > 0) 
and the trivial phase {5t < 0) can both be fully localized 
by disorder, because for either sign of St, the ground state 
|G) and top state |T) both belong to the same phase (we 
choose the convention that (2j —1, 2j) is a unit cell). This 
can be seen in the clean limit with At = 0. In momentum 
space H = (F{k)T^+(P{k)T^, and dis a nonzero 0(2) 

vector in the entire Id Brillouin zone with St y 0. For 
either sign of St, H and —H have the same topological 
winding number rii = ^ ^ dk FOkC^eab] thus |G) and 
|r) belong to the same phase. Based on our argument, 
all the finite energy states with either sign of St can be 
fully localized by random hopping At. The only states 
not fully localized in the two dimensional phase diagram 
tuned by e and St are located at the critical line St = 0. 
The critical line dt = 0 is in a “infinite-randomness” fixed 
point, and it can be understood through the strong dis¬ 
ordered real space renormalization group HMH]. 

Here we confirm the conclusions in Ref. [311110 that 
the finite energy density excited states of the Kitaev’s 
chain with St > 0 are still “topological”. Since the en¬ 
ergy level spacing between two eigenstates vanishes in 
the thermodynamic limit, the best way to determine if 
an excited state is topological or not is to compute its en¬ 
tanglement spectrum (the system is defined on a periodic 
Id lattice). And because the system is noninteracting, 
we will compute the single-particle entanglement spec¬ 
trum introduced in Ref. M for each excited state. The 
single-particle entanglement spectrum for the topologi¬ 
cal phase {St > 0) is shown in Fig.j^a), where two zero 


energy modes can be observed in the spectrum (corre¬ 
sponding to the Majorana zero modes at both entangle¬ 
ment cuts respectively). This topologically non-trivial 
feature persists for all energy eigenstates in the many- 
body spectrum, including the ground/top states and the 
finite energy density states in between. However at the 
critical line St = 0, as shown in Fig.j^b), the zero en¬ 
ergy modes are lifted by the long-range entanglement, 
and the single-particle entanglement levels become gap¬ 
less around = 0 which leads to the logarithmic scaling 
of the entanglement entropy. 

The Kitaev’s chain itself is just a free fermion model. 
But our argument indicates that under interaction, as 
long as |G) and |T) are still both in the topological su¬ 
perconductor phase, all of the excited states can still be 
localized. Such a generalization is justified given that 
the non-interacting Anderson localized states can be adi- 
abatically connected to the many-body localized states 
under interaction, as proven in Ref. [3]. 
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FIG. 1: Single-particle entanglement spectrum for many-body 
eigenstates of the random Kitaev’s chain, at (a) St = 0.5t and 
(b) St = 0. In both cases = 0.3t. We take a 128-site sys¬ 
tem with periodic boundary condition, which is partitioned 
into two 64-site subsystems for the entanglement calculation. 
e_B is the single-particle entanglement energy (s.t. the reduced 
density matrix pA = exp(—c^esc), as shown in Ref. 1191) . The 
spectrum of cb is shown as tanhe_B, and is calculated for sev¬ 
eral many-body eigenstates: including the ground and the 
top states and other 5 randomly picked finite energy density 
states, which are arranged in order of their energy density e. 
The shading denotes the standard deviation of the entangle¬ 
ment energy levels under a disorder average over the system. 
Of note are the topologically non-trivial, two-fold degenerate, 
zero energy modes throughout the entire spectrum e in the 
topological phase (a). 


2B. Modified Kitaev’s chain: critical states and 
interaction assisted localization 

In this subsection we consider a modified Kitaev’s 
chain: 

H ^^-{t- {-lyt'Fj + Atj) - hay (2) 
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where again Atj is random and t, t', h > 0. In this model 
tT| commutes with the Hamiltonian, which implies that 
any energy eigenstate will also be an eigenstate of cr|. 
In the clean limit, the ground state |G) of the system 
has (t| = I everywhere, and the fermions are in the triv¬ 
ial phase; in contrast, |T) must have cr| = — 1 every¬ 
where, and hence |T) is in the topological superconduc¬ 
tor phase. With disorder, both states can be localized, 
and their entanglement entropy shows the area-law scal¬ 
ing {i.e. S ^ const, for Id) as in Fig.j^a). But since 
the ground state and the top state belong to different 
SPT phases, based on our argument, there must be some 
finite energy density states which cannot be fully local¬ 
ized. In this model it is easy to visualize these delo¬ 
calized excited states. An excited state of the system 
has a static background configuration of cr| which does 
not satisfy = I. If we consider a random configura¬ 
tion of (t| that has the average = 0, then one can 
simply absorb ct| into the random numbers Atj, and 
the effective Hamiltonian for Majorana fermions 'jj reads 
HeS = ~ which is precisely the 

random hopping Majorana fermion model Eq. Q tuned 
to the critical point St = 0 . And according to Bef. IT^ ITKl 
the ground state of Hes (which is a highly excited state of 
the original Hamiltonian Eq. (U due to the h term) has a 
power-law correlation after disorder average, and its en¬ 
tanglement entropy scales logarithmically with the sub¬ 
system size: S ^ log i US]. So the delocalization happens 
right at the energy scale = —h J2j ^j = 0 . In deed 
our numerical calculation shows that as long as E^r ^ 0, 
the eigen states are all localized with area-law entangle¬ 
ment entropy as in Fig.[^a,b); but for E^ = 0, the eigen 
states are delocalized with logarithmically-scaled entan¬ 
glement entropy as in Fig.j^c). Thus the model Eq. ([^ 
cannot be fully many-body localized, which is consistent 
with our statement made in the introduction. 

The model Eq. has a time-reversal symmetry T : 
—>■ {—y^j and cr| —>■ cr|. It is known that with 
this time-reversal symmetry and without interactions, 
the Kitaev’s chain has Z classification [ 201122 ) : that is 
with an arbitrary number of flavors of Eq. |T) is al¬ 
ways a nontrivial topological superconductor, while |G) 
is always a trivial phase. However under certain flavor 
mixing four-fermion interaction [ 231124 ) . the classification 
of Kitaev’s chain with time-reversal symmetry reduces 
to Zg. Namely under this four-fermion interaction, for 
eight copies of Eq. |G) and \T) become the same triv¬ 
ial phase, which implies that there does not have to be 
any phase transition when increasing e, and all of the 
finite energy density excited states can be fully localized 
under the interplay between disorder and interaction. 

In model Eq. [^ the logarithmic entanglement entropy 
at the critical excited state comes from the long range 
effective hopping under renormalization group [iMni- 
We can assume that the four-fermion interaction on each 
site is random, then when and only when there are 8k 
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FIG. 2: Entanglement entropy Si vs log subsystem length 
log 2 i vs fermion energy E-, (energy of the first term in Eq. §) 
for various boson energies cr| = —hL, —hL/2,0 

(a,b,c) (second term in Eq. [^. Calculations are done on a 
random Majorana chain with L = 1024 sites, and the stan¬ 
dard deviation of Atj is CAt = t. States with Ea = 0 are 
the critical excited states which are delocalized. All states 
with different E-y at Ea = 0 have logarithmic entanglement 
entropy, and hence are delocalized. 


copies of Eq. [^ under interaction each site independently 
possesses a random set of many-body spectrum without 
degeneracy. Let SV be the typical energy level spacing 
of the interaction Hamiltonian on each site. To cre¬ 
ate entangled pairs between distant sites, the effective 
long-range coupling teff generated under RG must over¬ 
come the energy scale of SV to hybridize the many-body 
states. However the effective coupling strength actually 
falls rapidly with the distance [l3Hl5) as teff ~ te~'^, so 
the long-range coupling can only lead to exponentially 
small entanglement AS ^ {tes/SVy ~ {tjSVye~'^'^. 
Therefore even with weak interaction, all of the eigen¬ 
states are short-range entangled area-law states, and can 
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be fully localized. In contrast, without interaction, no 
matter what kind of fermion-bilinear perturbations we 
turn on in Eq. as long as these terms preserve the time- 
reversal symmetry defined above and the topological na¬ 
ture of jG) and |T), there must necessarily be some fi¬ 
nite energy density states which cannot be fully localized. 
Thus in this case interaction actually “assists” many- 
body localization, which is opposite from what is usually 
expected for weak interaction, in for example Ref. 1251 and 
is also different from the strong interaction reinforced lo¬ 
calization studied in Ref. [2SIB51 

Notice that this “interaction assisted localization” is 
only possible with 8k copies of the Kitaev’s chain with 
time-reversal symmetry. With 4 copies of the Kitaev’s 
chain, the spectrum on each site contains two sets of 
two-fold degenerate states even under interaction that 
preserves time-reversal, then the effective long-range cou¬ 
pling teS generated under RG will still lead to maximal 
entanglement between distant sites. A detailed RG anal¬ 
ysis about this will be given in another paper |29| . 

2C. Bosonic SPT states, Haldane phase 

Many bosonic SPT parent Hamiltonians can be written 
as a sum of mutually commuting local terms. For exam¬ 
ple, the “cluster model” for the Id SPT with Z 2 x Z 2 
symmetry m, the Levin-Gu model [30] and the GZX 
model [HI for the 2d SPT states with Z 2 symmetry, 
and the 3d bosonic SPT state with time-reversal sym¬ 
metry |32] are all a sum of commuting local operators; 
thus their ground states are a product of eigenstates of 
local operators m- SPT Hamiltonians written in this 
form are very similar to the “universal” Hamiltonian of 
MBL state proposed in Ref. |331 which is also a sum of 
mutually commuting local terms, because a MBL system 
has an infinite number of local conserved quantities. 

All of the idealized SPT models mentioned above have 
a Z 2 classification, and their ground and top states belong 
to the same SPT phase. Obviously there should be no 
phase transition while increasing energy density e. This 
statement is still valid with small perturbations which 
make these models nonintegrable as long as the nature of 
|G) and |r) are not affected by the perturbations. Thus 
these models (and their nonintegrable versions) can all 
be fully localized by disorder. 

However, some other bosonic SPT models can not be 
fully localized. In the following we will give one such 
example for the Haldane phase (SH jSi] : 

H = ■ ^^+1 + • • • (3) 
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Sj are spin-1/2 operators. The ellipsis includes perturba¬ 
tions that break the system’s symmetry down to a smaller 
symmetry (such as time-reversal or Z 2 x Z 2 ) that is suf¬ 
ficient to protect the Haldane phase, but do not lead 


to degeneracy in the bulk spectrum, namely only the 
boundary transforms nontrivially under symmetry. If 
the random coupling AJj is not strong enough to change 
the sign of J, then the ground state and top state of 
this model correspond to two opposite dimerization pat¬ 
terns of the spin-l/2s. Thus one of them is equivalent 
to the Haldane’s phase while the other is a trivial phase 
as long as we pick a convention of boundary. If we as¬ 
sume the random Heisenberg coupling A J is sufficient to 
localize most of the excited states, then there must be 
an unavoidable phase transition while increasing energy 
density e. According to our argument in the introduc¬ 
tion, this phase transition should behave just like an or¬ 
dinary quantum phase transition at zero temperature. It 
is known that the quantum phase transition between a 
Haldane phase and a trivial phase is a conformal field 
theory, and it is equivalent to a spin-1/2 chain without 
dimerization. With strong disorder, this quantum criti¬ 
cal point will be driven into the inhnite-randomness spin 
singlet phase mHni with a power-law decaying disorder 
averaged spin-spin correlation function and a logarithmic 
entanglement entropy m- 

2D. 2d interacting topological superconductor: 
critical states and interaction assisted localization 

In this subsection we will discuss the nonchiral 2d p±ip 
topological superconductor, i.e. p + ip pairing for spin-up 
fermions, and p — ip pairing for spin-down fermions. On 
a square lattice this TSC can be written in the Majorana 
fermion basis: 

H = ^ sin fca,-kr^cr"’sin 

k 

+ x‘_fer^(e - cos - cos ky)xk, (4) 

where = ±1 represents spin-up and down, while 
T'^ = ±1 represents the real and imaginary parts of 
the electron operator. Without any symmetry, this sys¬ 
tem is equivalent to the trivial state, i.e. its bound¬ 
ary can be gapped out without degeneracy. However, 
when 0 < e < 2, with a Z 2 symmetry which acts as 
Z 2 : X —)■ cr^x, the system is a nontrivial TSG. This sys¬ 
tem can also have another time-reversal symmetry, which 
is unimportant to our analysis. The boundary of this sys¬ 
tem reads: H = J dx X*(~*^a:^^)X; Z 2 : x ^ The 
Z 2 symmetry forbids any single particle backscattering at 
the boundary for arbitrary copies of the system, thus the 
p ±ip TSG with the Z 2 symmetry has a Z classification 
without interaction. 

Without any interaction, for n—copies of the p A ip 
TSG, |G) and |T) belong to different SPT phases. This 
is because for either spin-up or down fermions, the Chern 
number of |G) and |T) are opposite. And because the sys¬ 
tem has a Z classification, |G) and |T) must belong to 
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different SPT states. Using our argument in the intro¬ 
duction, this implies that under disorder that preserves 
the Zi symmetry, there must be some finite energy den¬ 
sity states which cannot be fully localized. This is not 
surprising, considering that even at the single particle 
level there are likely extended single particle states un¬ 
der disorder. The existence of extended single particle 
states is well-known in integer quantum Hall state [36], 
and recently generalized to quantum spin Hall insulator 
with a Z 2 index |S71|35|. 

The situation will be very different with interactions. 
Once again a well-designed interaction will reduce the 
classification of this p ± ip TSC from Z to Zg . 

Namely n—copies of Eq. [^ is topologically equivalent to 
(n-I-8fc)—copies. This implies that under interaction |G) 
and \T) actually belong to the same phase when n = Ak. 
Thus when n = 4fc, the phase transition in the noninter¬ 
acting limit will be circumvented by interaction above a 
certain critical value. Thus once again interaction assists 
MBL in this case. When n = 8, |G) and |T) are both 
trivialized by interaction, namely interaction can adia- 
batically connect both states to a direct product of local 
states. When n = 4, Ref. H51 showed that interaction can 
confine the fermionic degrees of freedom, and drive four 
copies of the p ± ip TSC into a 2d bosonic SPT state 
with Z 2 symmetry, which as we discussed in the previous 
section, can also be fully many-body localized. 

Please note that in the noninteracting limit the quan¬ 
tum phase transition between 2d TSC and trivial state 
is described by gapless (2 -|- l)d Majorana fermions, and 
since a weak short range four-fermion interaction is ir¬ 
relevant for gapless (2 -|- l)d Dirac/Majorana fermions, 
only strong enough interaction can gap out the quantum 
phase transition. Thus unlike the Id analogue discussed 
in section2B, we expect that in this 2d system only strong 
enough interaction can “assist” disorder and localize all 
the excited states even for n = 4fc. 


3. SUMMARY 

In this work we propose a simple rule to determine 
whether a local Hamiltonian with symmetry can be 
many-body localized. Since MBL is a phenomenon for 
the entire spectrum, we need to start with a lattice 
Hamiltonian for our analysis. Therefore the low energy 
field theory descriptions and classification of SPT states 
such as the Chern-Simons field theory [44] and the non¬ 
linear sigma model field theory will not be able to 
address this question. Instead, our argument is based 
on the nature of the ground and top states of the same 
lattice Hamiltonian. Our argument is general enough, 
that it can be applied to both free and interacting sys¬ 
tems, bosonic and fermionic SPT systems. And coun¬ 
terintuitively, we found that because interactions change 
the classification of fermionic topological insulators and 


topological superconductors, in some cases interactions 
actually assists localization, rather than delocalization. 
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